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Abstract 

Recent LHC results concerning full jet-quenching in Pb Pb collisions have been 
presented in terms of a jet asymmetry parameter, measuring the imbalance between 
the transverse momenta of leading and subleading jets. We examine the potential 
sensitivity of this distribution to fluctuations from the heavy-ion background. Our 
results suggest that a detailed estimate of the true fluctuations would be of benefit 
in extracting quantitative information about jet quenching. We also find that the 
apparent impact of fluctuations on the jet asymmetry distribution can depend sig- 
nificantly on the choice of \ow-pt threshold used for the simulation of the hard pp 
events. 

In the quest to understand the properties of the medium generated in high-energy 
heavy-ion collisions, the past decade has seen extensive study of medium-induced modifi- 
cations to the production of high transverse momentum objects [1]. It has been conclu- 
sively established at RHIC that the spectra of high-momentum hadrons are significantly 
suppressed, by a factor of Raa — 0.2 relative to the appropriate rescaling of the pp spectra. 
This effect is generally attributed to their (or their originating parton's) interaction with 
the medium. 

Recently, significant attention has been directed to jets. Compared to hadrons, jets are 
interesting because, at least in pp collisions, there is a closer, and perturbatively quantifi- 
able, connection between a jet's momentum and that of its originating parton. STAR [2] 
and PHENIX [3] have presented first (preliminary) measurements of full jets produced in 
AuAu collisions with transverse momenta in the 20 — 40 GeV range and found that jet 
spectra are also suppressed, though by a potentially more modest factor than for hadrons. 
Recently, ATLAS has published studies of the correlations between the momenta of the 
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Figure 1: A simulated pp event from Pythia 6.423 (centre-of-mass energy ^/s = 2.76 TeV; 
the missing transverse energy was zero). We find that for 1 in every 300 events with a jet 
with pti > 100 GeV, the second hardest jet has pt2 < Pti/^- A more accurate estimation 
of this number would benefit from the combination of 2-jet, 3-jet, 4-jet, . . . samples, using 
for example one of the multijet matching methods reviewed in [8]. 



two leading jets, with the striking observation that a significant fraction of events show 
a strong imbalance between the pts of the leading jet and the first subleading jet on the 
opposite side of the event. CMS has shown similar preliminary results in Ref. [Sjj^and first 
phenomenological interpretations have been given in Ref. [7] . 

Dijet imbalances can occur also in normal pp events, due to emission of multiple gluons 
(cf. the simulated Pythia event shown in fig. [T]), but they are quite rare. To quantify how 
much more often they arise in Pb Pb collisions, ATLAS and CMS have shown distributions 
of the jet asymmetry Aj 

Eti — Et2 X 

expressed in terms of the transverse energies of the leading and subleading jets, respectively 
E-Ti and Et2- The main quantitative evidence for jet quenching comes in the form of a 
significant enhancement of the asymmetry in the region around Aj ~ 0.4 (fig. 3 of [4j and 
p. 26 of [5]). 

In extracting the distribution of Aj, the experiments must contend with the fact that 
each jet may be contaminated with O (100 — 150 GeV) of transverse momentum from the 
medium particles, usually referred to as the background]^ To calculate the Aj for a given 



^Subsequently published in [6]. 

^The distinction between medium particles and jet particles is not necessarily very legitimate physically, 
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dijet event, each jet's momentum is corrected for the expected level of background activity 
in the jet, usually estimated from the activity elsewhere in the event, preferably at similar 
rapidities (see e.g. jUIIUlll]). Such a correction cannot, however, account for the fact that 
the background fluctuates from point to point within the event (even at the same rapidity), 
so that the momentum subtracted from the jet will inevitably differ from the background 
actually present in the jet; nor does it account for fluctuations in the detector's response 
to the background and jet particles. 

Fluctuations are of course a common issue for jet measurements even in pp collisions, 
notably due to randomness in the response of calorimeters. However two novelties may 
be relevant concerning fluctuations for heavy-ion collisions. Firstly the LHC heavy-ion 
medium has only just been produced and it is probably fair to assume that its fluctuation^ 
are less well understood than those of the detectors, which have been the object of study for 
many years. Secondly, the absolute size (i.e. in GeV) of detector fluctuations scales roughly 
as the square-root of the jet energy, meaning that they are less important for low-pt jets 
than for high-pj jets, whereas background fluctuations are probably largely independent of 
the jet's pt- 

This last point is relevant because of the way in which fluctuations can affect the 
Aj distribution. The experimental analyses of the Aj distribution select events in which 
the leading jet passes some high-p^ cut, say pt > 100 GeV. Events with a genuine high- 
Pt jet are rare. There are many more low-pf dijet events and in some small fraction of 
cases the background under one of the jets may fluctuate upwards causing the jet to 
pass the high-pt cut. Such events will naturally have a large jet asymmetry, since there 
is no reason for the balancing jet to also have a positive background fluctuation. The 
relative contributions of different classes of events depends on the interplay between the 
rareness of large background fluctuations and the rareness of high-pj jet production as 
compared to low-pt jet production. While one can in principle estimate the potential 
severity of this problem from Monte Carlo simulations, it is debatable whether reliable 
enough descriptions of the PbPb medium produced at high energy exist. Guidance from 
experimental measurements is therefore paramount. 

One parameter that is indicative of the size of the fluctuations in the reconstructed jet 
Pt is their standard deviation, which we call cxjet- ATLAS [H] has presented preliminary 
results for the fluctuations from one calorimeter tower to the next. If scaled up by the 
square-root of the number of towers in a jet (about 50 for an R = 0.4 jet with towers of 
size 0.1 X 0.1 in rapidity and azimuth), it would suggest a value cijet — 8.5 GeV for the 
most central set of events. On the other hand, the scaling of the tower fluctuations by 
the square-root of the number of towers may not be a safe way of extrapolating tower 
fluctuations to a jet, insofar as the background could well have local correlations (there is 
no clear reason for the correlation length of such fluctuations to necessarily be smaller than 
the calorimeter tower size). Furthermore there can be other factors that contribute to a 

however it may still make sense to think of an expected level of background transverse momentum. 

■^Including their standard deviation, correlations from point to point within the event, non-Gaussianities, 
etc. 
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degradation of resolution, such as back reaction [12] and fluctuations in the event-by-event 
(or calorimeter-strip) estimation of the background level (as discussed in sections 3.5 and 
A.l of [ID]). 

Another way in which one may attempt to deduce the level of the fluctuations is from 
a preliminary inclusive jet spectrum for — 100% centrality (p. 41 of [H]), which displays 
a region of near Gaussian pj-dependence for pt < 50 that is strongly suggestive of an 
origin due to fluctuations, and compatible with ajet — 14 GeV. One would then expect the 
corresponding ajet for — 10% central events to be somewhat larger. 

Neither of the above estimates, made in the original version of this article, is ideal. One 
criticism that has been made is that the tower fluctuations do not include full calibration. 
As for the inclusive jet spectrum, it mixes a genuine jet spectrum in with the fluctuations. 
However, subsequently, support for crjct values > 15 GeV has appeared in the form of 
a preliminary measurement of the background fluctuations for charged-track jets from 



ALICE [13]. It is discussed in detail in Appendix A. 2 

To provide simple insight into the impact on the dijet asymmetry from various values 
of (Tjet, we have carried out the following "toy" analysis. We generate jet events with 
Pythia [H] (version 6.423, DW underlying-event tune p!5|). To mimic the effect of residual 
fluctuations following background subtraction, we then add to the pt of each jet a random 
fluctuation, generated according to a Gaussian distribution with mean and standard 



deviation cxjet, independently of the jet pt (details are given in appendix A.l). These 
choices correspond to a perfect estimate of the average background that needs subtracting 
in each event, with crjet encoding the combination of all sources of fluctuations, whether 
intrinsic fluctuations of the background, or fluctuations in the detector's response to it. We 
select events in which the leading jet has pt > 100 GeV, the subleading jet has pt > 25 GeV, 
both have rapidities \y\ < 2.8 and are separated in azimuth by |A0| > 7r/2 and for these 
events plot the corresponding distribution of Aj, similarly to the ATLAS analysis [1]. 

The filled black points in fig. |2]show our results for four different values of crjct. One sees 
a clear distortion of the Aj distribution as cXjet is increased, reminiscent of the pattern seen 
by ATLAS and CMS with increasing centrality. One key element of our simulation is that 
in generating the filled black points we chose a fairly low minimum pt cut, p™™ = 30 GeV, 
for the underlying Pythia 2 — 2 scattering, and also verified that further lowering this 
cut made no difference for our values of cijet. With a larger choice, pf = 70 GeV (shaded 
region) which would be perfectly adequate for low crjet, one notices that a significant part 
of the effect of the background fluctuations can be missed for larger crjet- This leads to the 
obvious implication that the choice of can play an important role, especially if cXjet 
happens to be large (or, as we have also found, if there are significant non-Gaussianities in 
the fluctuation^. 

A complementary investigation into the impact of fluctuations can be obtained by 
embedding Pythia events into a simulated PbPb background. A similar investigation was 
carried out by ATLAS, embedding events into PbPb events as simulated by an ATLAS- 



''We understand that this was the value used in refs. [H [TT] 
""Significant non-Gaussianities have been observed in [16]. 
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Pythia with Gaussian smearing 
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Figure 2: Simulated distribution of Aj and A<^, as obtained when smearing the pt of jets 
from Pythia 6.4 (DW tune [I5]) by an amount cxjet- None of the results in this figure involved 
jet quenching. Four different ajet values are shown, and for each plot there are results 
from Pythia simulations with two different generation cutoffs on the 2 — )■ 2 scattering, 
pmm _ 20 pmin _ -jrg QgY^ gQ g^g illustrate its impact. The results labelled "pp" 

reference always correspond to p™™ = 30 GeV with no smearing. Jet clustering has been 
performed with the anti-fct algorithm ^17j with R = 0.4, as implemented in Fast Jet |18j . 



specific version of HIJING \i9\. Our analysis will differ in that we study HYDJET [20] 
rather than HIJING and use also a lower pt cutoff for the Pythia events. The tune we 
use for HYDJET]^ gives an average background level of 210 GeV per unit area for — 10% 
centrality and \ri\ < 2.8, compatible with the average jet contamination found by ATLAS, 
and an average charged-particle multiplicity for — 10% centrality of 1400 for Ir^l < 0.5, 
which is reasonably consistent with that measured by ALICE [21j (further comparisons are 



discussed in appendix A. 2). HYDJET's simulation of quenching has been turned off, to 
avoid the potential confusion that might arise from the quenching of hard jets associated 
with the PbPb simulation rather than with the embedded Pythia event (quenching has 
only a modest 5 — 10% effect on the HYDJET fluctuations). Since detectors can have an 
impact on fluctuations, we have also processed the events through a simplified calorimeter 



^The tuning parameters used to simulate LHC events at -y/s — 2.76 TeV with HYDJET vl.6 have been 
extrapolated between the 200 GeV (RHIC) and 5.5 TeV (LHC at designed energy) values used in [TU] 
(footnote 7), namely nh = 25600, ylfl = 3.9, ytfl = 1.46 and ptmin = 7.54 GeV. Quenching effects have 
been switched off by setting nhsel = 1. The embedded events come from Pythia version 6.423, tune DW, 
run at y/s = 2.76 TeV. 
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simulationj^ To subtract the background from jets we have taken the area/median method 
of [121 122] , using, for the estimation of the background density, a (rapidity) Strip Range of 
half- width 0.8 centred on the jet to be subtracted, as described in more detail in [TU]. This 
method should perform similarly to the ATLAS method of background subtraction. With 
this setup, for collisions in the 0-10% centrality range, we find fluctuations per unit area 
of about 23 GeV corresponding, for anti-fct jets of radius R = 0.4, to an expected crjet of 
16 GeV and a measured crjct of 17 GeV. 

The results we obtain from the HYDJET-|-Pythia simulations are presented in Fig. [3] 
for four centrality ranges. The empty circles labelled "pp" reference correspond to plain 
Pythia results as for Fig. |2} The filled black points and shaded histogram correspond to 
our embedding in HYDJET events and differ only by the of the underlying Pythia 
2 — 2 scattering: 10 GeV has been used for the former and 70 GeV for the latter]^ 

The evolution of the Aj distribution with increasing centrality in HYDJET displays 
a pattern similar to that observed for the Gaussian smearing with increasing cjjet- If 
anything, the distortion of the Aj distribution for — 10% central HYDJET collisions is 
slightly more pronounced at large Aj than with the highest Gaussian smearing we used, 
despite the smaller cxjet value from HYDJET. This could perhaps be a consequence of non- 
Gaussianities in its fluctuations. The HYDJET results also confirm the importance of the 
choice of p™™ cut on the 2 — )■ 2 scatters. 

While the above results suggest that fluctuations could be of relevance in interpreting 
the Aj distributions, one should not forget that the experiments have studied observables 
intended to signal the possible presence of important effects from fluctuations. One such 
observable is the fraction of energy inside a core of i? = 0.2 within the jet. A fluctuation 
that enhances the leading jet's pt would not necessarily be close to the centre of the jet 
and so should on average reduce the core energy fraction]^ Preliminary data from ATLAS 
(p. 34 of [H]) show a stronger reduction in core energy fraction with increasing centrality 
than in the ATLAS HIJING simulations. In our HYDJET simulations, the core energy 
fraction decreases yet more rapidly, which at first sight suggests that its fluctuations could 
be excessive. On the other hand, we find that the agreement in absolute value is better for 
central collisions than for peripheral collisions, complicating the interpretation]^ Another 
cross-check on fluctuations comes from the Aj distribution for jets with R = 0.6 (e.g. p. 48 



^Charged particles with pt < 0.5 GeV are first removed, and the remaining particles are put on a 
calorimeter of size 0.1 x 0.1 extending up to I77I =5 with uncorrelated Gaussian fluctuations of standard- 
deviation 0.8/y/E in each tower and a GeV tower threshold. This simple procedure gives a resolution that 
is 5 — 10% better than the true ATIAS calorimeter resolution; the comparison and other relevant points 
are discussed in appendix A. 3 The number quoted above for the average energy flow and fluctuations are 



those obtained at calorimeter level. 

^HYDJET itself generates many additional pp 2 2 scatterings for each heavy-ion collision, each with 
ptmin = 7.54 GeV. When embedding a jet event with a 10 GeV cutoff, in most cases the two hardest jets 
actually originate from these additional HYDJET pp scatterings. 

^This though is not entirely trivial, because the fluctuation may itself displace the centre of the jet. 
Furthermore any quenching of the leading jet may also reduce the core energy fraction. 

-•^^For the subleading jet, the centrality dependence is very similar for data and HYDJET, but the data 
are systematically about 0.15 below HYDJET. 
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Pythia embedded in HYDJET 
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Figure 3: Simulated distribution of Aj and A<^, as obtained when embedding Pythia 
events in a PbPb background described by HYDJET 1.6. None of the results in this figure 
involved jet quenching and the results obtained with HYDJET include a simple calorimeter 
simulation. Four different centrality regions are shown as indicated in the plots on the top 
row. For each plot there are results from Pythia simulations with two different generation 
cutoffs on the 2 — )■ 2 scattering, p™™ = 10 GeV and p™™ = 70 GeV, so as to illustrate 
its impact. The results labelled "pp" reference always correspond to those of Fig. [2j Jet 
clustering has been performed with the anti-fct algorithm [T7] with R = 0.4, as implemented 
in Fast Jet [18] and the heavy-ion background subtraction has been performed as described 
in [To] with the background density estimated using a Strip Range of half- width 0.8 centred 
on the jet being subtracted. 

of [H]). Since fluctuations should increase for a larger R, one would expect them to lead 
to an enhancement of the high Aj part of the R = 0.6 distribution. Vacuum QCD (and 
jet quenching) are expected to act in the opposite direction. The (unquenched) HYDJET 
simulation shows a fairly complicated behaviour however: the large Aj (> 0.4) part of 
the distribution barely changes in going from R = 0.4 to R = 0.6, while the distribution 
increases for Aj = 0.2 (and decreases for Aj near zero). In contrast, the preliminary data 
decrease at large Aj and, within the (large) errors, barely change for moderate and small 
Aj, suggesting, possibly, some non-trivial interplay between an effect such as quenching 
and fiuctuationsPI 

^^Data from CMS [5] on momentum flow in tracks, which appeared subsequent to the first version of 
this article, also indicate some genuine component of quenching. 
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If fluctuations are relevant, as they seem to be, then it is probably advantageous to 
attempt to unfold their effect, and/or to reduce their impact by raising the jet pt thresholds. 
Additionally it may be of interest to investigate methods to suppress fluctuations (the 
method of [9j as used by CMS [B], or filtering/trimming/pruning [231 EH [2^). Nevertheless 
one should be aware that such methods introduce potential biases of their own, as has been 
found in earlier studies [21 [10], and it then becomes important to quantify the interplay 
between any quenching and the noise reduction method,^ 

To conclude, we have found that fluctuations can significantly affect the dijet asymmetry 
A J measured in [4]. A precise estimate of the contribution of fluctuations is therefore 
important to be able to quantify the degree of quenching that is present in the data. A first 
direct estimate of these fluctuations has appeared in preliminary form [13j since the original 
version of our article. It shows fluctuations that are well reproduced by our HYDJET 
simulations and consistent also, therefore, with the upper end of the range that we explored 
in the toy model. Quantitative use of the Aj data to constrain quenching therefore probably 
requires that the potential bias due to fluctuations (or from any techniques used to suppress 
them) be accounted for. 
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A Appendix 

Since the appearance of the first version of this article, a number of questions have been 
raised concerning our analysis. Moreover new experimental results have become available. 
We address these questions here, also in light of the new data. 



A.l Details of our Gaussian smearing procedure 

In the discussions that followed the appearance of our article, one question that arose 
concerned the details of our Gaussian smearing procedure, which we had originally omitted 



Such methods discard low-momentum components of the jets, exploiting the fact that the background 
is almost entirely made of low-momentum particles, while for a pp jet only a small fraction of its total 
momentum is contained in low-momentum components. In the presence of quenching, however, a larger, 
but unknown, fraction of the jet's energy may be concentrated in low-momentum components. Discarding 
these components is then not without risks. Special care should also be taken with infrared or coUinear 
unsafe seeded jet algorithms, as quenching may cause a jet's high central energy density (the seed) to be 
redistributed over a broader region of the calorimeter. Further concerns, specific to the method used by 



CMS, are discussed in appendix A. 4 
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for brevity. The procedure is the following: we generate Pythia events (pp, y/s = 2.76 TeV, 
DW tune); to these we add a low density of "ghost" particles (10 per unit area), which 
have infinitesimal momentum and serve to ensure that there are no substantial empty 
regions in the event; we then run the anti-fct algorithm with i? = 0.4 on the combination of 
Pythia and ghost particles. For each jet that is found (however small its momentum, and 
independently of whether it involves Pythia particles or just ghost-particles) we add to it 
a random pt, chosen according to a Gaussian distribution of zero mean and dispersion cxjet- 
We then consider only the two hardest jets, requiring that they both be above the 25 GeV 
threshold. 

A criticism was made [26] that only "hard" jets should be fluctuated and that if this 
choice was made then the deformation of the Aj distribution would differ substantially 
from that seen in fig. [2j The criticism may have its origins in the fact that one of the LHC 
experiments, CMS, uses a seeded jet algorithm. Depending on the seed threshold, then it 



might indeed be the case that a hard core is a prerequisite for a jet to be found, ^ For a 
non-seeded algorithm such as anti-Zc^, fluctuations anywhere in the event may cause a jet 
to be found, even if there was no corresponding genuine underlying hard jet. This has to 
be accounted for in the choice of original jets that receive fluctuations. We note that had 
we not fluctuated zero- and low-pt jets we would not have reproduced the broadening of 
the low A0 tail in fig. [2] (which for cXjet = 20 GeV is comparable to ATLAS's results for 
- 10% centrality). 

One criticism that would be legitimate is that some fraction of the anti-kt jets have a 
very small area and should therefore be subject to smaller fluctuations, whereas we assigned 
identical fluctuations to all jets. We have examined the impact of a modified version of 
our procedure, in which the fluctuations scale in proportion to the square-root of the jet 
area. The results show slightly larger fluctuation-induced asymmetries, perhaps because a 
modest fraction of anti-/ct jets have areas larger than nR^. 



A. 2 Fluctuations in ALICE data and HYDJET simulations 

To date, the most direct experimental constraint on cXjet comes from an ALICE result 
presented in a talk [13] a few weeks after the appearance of our note. This involves a 
measurement of the resolution in the reconstruction of jets from charged tracks only, tested 
via the embedding of single-track jets. The background subtraction method used there is 
the area/median method of Fast Jet. We can directly compare to it by applying our analysis 
procedure on the charged tracks from HYDJET, with a single additional embedded hard 
track. 

That comparison is given in fig. |4], which shows the difference between the reconstructed 
Pt of the jet that contains the track and the pt of the track itself. The agreement between 

^■^CMS actually has a quite low seed threshold, of 1 GcV 0, and so in practice the heavy-ion background 
may produce a high multiplicity of seeds in its own right. The choice of seed threshold is delicate precisely 
because too low a value brings in many "fake" seeds and too high a value may cause genuine jets to be 
missed, especially if quenching changes the relative proportion of high-pt particles in the jet. Additionally, 
seeds are intrisincally coUinear unsafe. 
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Figure 4: Comparison of our HYDJET simulation to ALICE data [13], showing the dis- 
tribution of Api = Pi,rec — Pt,track i-G. the difference in transverse momentum between a 
reconstructed charged-track jet (anti-fcj, R = 0.4) and the single embedded charged track 
contained within it. Apt is calculated after background subtraction with the Fastjet me- 
dian/area method. The jets have been reconstructed using charged particles with |?7| < 0.8 
and pt > 150 MeV, whose masses have been set to zero by rescaling their energy. The 
background estimation in the HYDJET case used the kt algorithm [27| with R = 0.4, with 
particles and ghosts (of area 0.01) up to \y\ = 0.8 and a global range up to \y\ = 0.4 (ex- 
cluding the two hardest jets within \y\ < 0.8). We understand that these choices coincide 
with those made by ALICE. 



10 



the ALICE and HYD JET results is striking, especially considering that our HYD JET simu- 
lation had not been directly tuned to LHC data. Quantitatively this can seen by comparing 
the crjet(chg) values as obtained from the ALICE data with those from the simulation. We 
find crjet(chg) ^ 11-5 GeV from the ALICE data and 11.4 GeV in the HYDJET simulation]^ 
Several factors intervene in comparing these results to the crjet value of 17 GeV discussed 
in the main text for full jets embedded in HYDJET. Firstly, in going from the embedding 
of a single track to that of a full jet, an extra increase of a few percent may arise (e.g. due 
to back-reaction [10]). More importantly, assuming that a fraction f^hg — 0.6 of the energy 
flow is carried by charged particles, one may expect an additional factor that is somewhere 
between ~ 1/ \J fchg — 1-3 (taking charged and neutral components to be uncorrelated) and 
~ l//chg — 1.7 (for 100% correlated charged and neutral components). Finally there will 
be additional fluctuations from the calorimetric nature of jet measurements. Combining 
all these factors gives a result that is consistent with the cijet — 17 GeV quoted in the main 
text. 

Ultimately therefore, the ALICE results support the choices that we made in estimating 
the possible order of magnitude of fluctuation-induced contributions to the jet-asymmetry 
measurement. Note however that other issues also need to be taken in to account at the 
1 — 2 GeV level (for example detector effects). An analysis at this level of accuracy can 
therefore probably only be performed in conjunction with a fu 11 detector simulation. 

A. 3 Suitability of calorimeter simulation 

One objection that has been raised concerning our HYDJET results is that our toy 
calorimeter simulation suffered from an overly pessimistic resolution. To obtain a detailed 
description of jet response, ideally many effects need to be taken into account, including 
tower thresholds, different kinds of noise term that are independence of energy, scale as 
l/\fE and as 1/-E, different responses to photons and charged hadrons, the effect of mag- 
netic fields, and the degradations of resolution due to detector elements that lie between 
the interaction point and the calorimeter. A full simulation of all these effects goes beyond 
the scope of this article, therefore we made a simple approximation involving towers of 
size 0.1 X 0.1 with a resolution of 0.8/-\/E, together with the removal of charged tracks 
with pt < 0.5 GeV (on the grounds that they would be bent away from the calorimeter 
by the magnetic field). The 0.8/ y/E term is larger than the corresponding term for the 
ATLAS hadronic calorimeter (~ O.S/a/E), however this difference is justified by the fact 
that we neglect many other sources of detector fluctuation. The validation of this state- 

There is a small residual shift between the ALICE and HYDJET results. When it is taken into 
account, the visual agreement in fig. |4] becomes perfect. The precise origin of this small difference has not 
been identified, but small differences in the subtraction procedure can lead to such shifts without affecting 
the crjet(chg) value (TU]. Note that the crjct(chg) can, however, be affected by the choice of ghost area. A 
smaller area can reduce it by 0.5 GeV. We have verified that fig. |3] remains essentially unchanged when 
using a smaller ghost area. Note that the crjot(chg) values quoted by ALICE were obtained from Gaussian 
fits to the Apt < region. Such fits give crjot(chg) values that are about 2 GeV smaller than results derived 
from the full distribution, which are more relevant to our studies here. 
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Figure 5: Comparison of jet resolution obtained with our toy calorimeter to the measured 
ATLAS jet resolution 



ment is given through fig. |5} which compares the jet resolution that we have observed with 
our calorimeter (matching calorimeter jets within AR < 0.25 of the two hardest particle- 
level jets in an event) to that found by ATLAS. One sees that rather than overestimating 
calorimeter fiuctuations, we actually underestimate them by 5 — 10% over the full range of 
Pt in which jet resolution measurements exist. Note also that ATLAS has shown that its 
calorimeter response is similar for jets and heavy-ion backgrounds (p. 37 of [llj). 

A further potential issue comes from the fact that an experiment's magnetic field, 
by bending particles differently according their pt, may smear any angular correlations 
that are present in the background fluctuations. Our calorimeter simulation did not by 
default consider the impact of a magnetic fleld (other than by removing charged particles 
with Pt < 0.5 GeV). We have, however, verifled that if one accounts for the effect of 
a magnetic fleld on the charged-particles' azimuth when they reach the calorimeter, the 
Fast Jet median/area estimate of fluctuations is reduced by an amount of order 1%. This 
value was obtained with a conflguration similar to that used by ATLAS, i.e. a longitudinal 
magnetic fleld of strength B = 2 T, with a calorimeter at a radius of 1.5 m from the beam. 
This result gives us confldence that our original approximation was not unreasonable. 

If is to be noted that magnetic flelds could also have an impact on measurements of 
quenched jets (independently of background- related issues), insofar as quenching may alter 
the relative fractions of different momentum components within a jet and therefore also 
give an appearance of additional jet broadening. 
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A. 4 The CMS jet-reconstruction procedure 

It was not the purpose of this article to discuss the CMS results [U], insofar as they appeared 
subsequent to the first version of this article. However, certain criticisms of our work have 
been made based on results presented by CMS, and so deserve comment. 

An objection that has been raised is that CMS explicitly show jet resolution results that 
are incompatible with the values we have investigated. For example, fig. 4f of [6J indicates 
that 40 GeV jets (0 — 10% centrahty) have about 20% resolution, i.e. a ajct — 8 GeV, 
well below the largest values that we have discussed here. In this context, however, it 
is important to be aware that CMS's jet reconstruction procedure differs substantially 
from atlas's. One of the differences is that the method used to subtract the heavy- 
ion background involves a noise reduction technique [9] that estimates the mean /i and 
standard deviation Utower of the calorimeter towers' pt deposits (after excluding hard jets) 
and then subtracts /i + atower from each tower's pt- Only towers that are positive are then 
retained. 

To understand the potential benefits and biases of this method, let us make the as- 
sumption that /i and cxtower are well determined (there are a number of complications in 
their practical determination, e.g. with respect to the exclusion of hard jets, which may 
degrade the performance of the method with respect to the analysis that follows). Let us 
also assume purely Gaussian noise (again, this is probably optimistic). Then the fraction 
of towers retained is 

1 _ 

dx^= e ^^- 0.1587, (2) 



where x = pt^tower — 1^- Since all of these towers are positive, they induce a systematic offset 
in the overall reconstructed jet momentum 

(t — rr, ) - 

(^pSD = A^towcr(<5pSwer) = ^^tower / dx^-j=^^e ^ 0.0833 CTtowcr A^tower , 

(3) 

where Atower is the total number of towers that are contained in a jet (vri?^ divided by the 
tower area, which is 0.087 x 0.087, i.e. A'towcr ^ 104 with R = 0.5 as used by CMS). For 
Ctower — 1 — 2 GeV, this would correspond to an 8 — 16 GeV bias. 

The residual fluctuations after the noise suppression can be estimated as 



[CF: 



no.e-suppressed,2 _ [((5p— J^) - (4a) 



jet / 



tower 



7 V-^ tower; WZ2 I j- noise \2 

c?2;^== e ^^tower - {Spt^^^^J 

Ttowcr V ■^^''"tower 

2 



(4b) 

(0.262 (Ttower)'Artower, (4c) 



i.e. a 74% reduction in the noise as compared to subtraction procedures without noise 
suppression, which just give af^^ ~ cr^o^gj, Atower- This is a potential strength of the CMS 
noise-reduction technique, but it comes at the price of introducing a significant pt offset, 
eq. (|3|. Since noise-reduction has such a large impact on the fluctuations it is not possible 



13 



to draw any conclusion about ATLAS results based on the quoted CMS jet resolution after 
noise-suppression. 

In practice, in simulations that embed hard pp jets in a heavy-background, the offset 
of eq. ^ will not be directly seen. This is because embedded hard jets are not just pure 
noise, but involve some number of towers that are far above the + atower threshold. These 
towers will all be retained. However, relative to the original pp towers, there will now be 
an offset of — atowcr for each of these towers. Defining the the calorimeter "occupancy" of 
a normal pp jet to be / (i.e. a pp jet contains on average /A^tower active towers), then the 
total offset from this effect will be 

(5pS)--/A^toweratower. (5) 

In the the limit in which / is small its impact can be neglected in eqs. 
QCD vacuum jets we find that it has a value / ~ 0.1,^ 

Taking into account this effect and the offset of eq. ([3]), gives us an overall offset of 

(5p?]:r") = (^pSD + ('^pS) ^ (0.0833 - /)Ar,„_atower • (6) 

Thus the net bias is small, but only due to a fortuitous cancellation between two effects 
with very different physical originsj^ Both of these effects are at the same 10% level as 
the overall impact of quenching (e.g. fig. 12 of p]). Furthermore, it is not immediately 
clear that / for quenched jets is the same as for pp jets and one may even expect it to 
be larger, especially for the away-side jet, leading to an additional negative offset for that 
jet, thereby artificially enhancing the asymmetry. Therefore any quantitative analysis of 
quenching in heavy-ion collisions that relies on noise reduction should also perform an 
analysis of systematic errors due to any biases associated with potential medium-induced 
modifications of /. 

The discussion also implies that it is difficult to relate the large number of supporting 
plots by CMS (e.g. track-properties of calorimeter jets) directly to the ATLAS case. 



3|4). For hard 
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